Efficient and reproducible generation of human iPSC-derived cardiomyocytes using a stirred bioreactor.

In the last decade human iPSC-derived cardiomyocytes (hiPSC-CMs) proved to be valuable for cardiac disease modeling and cardiac regeneration, yet challenges with scale, quality, inter-batch consistency, and cryopreservation remain, reducing experimental reproducibility and limiting clinical translation. Here, we report a robust cardiac differentiation protocol that uses Wnt modulation and a stirred suspension bioreactor to produce on average 124 million hiPSC-CMs with >90% purity using a variety of hiPSC lines (19 differentiations; 10 iPSC lines). After controlled freeze and thaw, bioreactor-derived CMs (bCMs) showed high viability (>90%), interbatch reproducibility in cellular morphology, function, drug response and ventricular identity, which was further supported by single cell transcriptomes. bCMs on microcontact printed substrates revealed a higher degree of sarcomere maturation and viability during long-term culture compared to monolayer-derived CMs (mCMs). Moreover, functional investigation of bCMs in 3D engineered heart tissues showed earlier and stronger force production during long-term culture, and robust pacing capture up to 4 Hz when compared to mCMs. bCMs derived from this differentiation protocol will expand the applications of hiPSC-CMs by providing a reproducible, scalable, and resource efficient method to generate cardiac cells with well-characterized structural and functional properties superior to standard mCMs.

Functional dissection of human cardiac enhancers and noncoding de novo variants in congenital heart disease.

Rare coding mutations cause ∼45% of congenital heart disease (CHD). Noncoding mutations that perturb cis-regulatory elements (CREs) likely contribute to the remaining cases, but their identification has been problematic. Using a lentiviral massively parallel reporter assay (lentiMPRA) in human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), we functionally evaluated 6,590 noncoding de novo variants (ncDNVs) prioritized from the whole-genome sequencing of 750 CHD trios. A total of 403 ncDNVs substantially affected cardiac CRE activity. A majority increased enhancer activity, often at regions with undetectable reference sequence activity. Of ten DNVs tested by introduction into their native genomic context, four altered the expression of neighboring genes and iPSC-CM transcriptional state. To prioritize future DNVs for functional testing, we used the MPRA data to develop a regression model, EpiCard. Analysis of an independent CHD cohort by EpiCard found enrichment of DNVs. Together, we developed a scalable system to measure the effect of ncDNVs on CRE activity and deployed it to systematically assess the contribution of ncDNVs to CHD.

Molecular and Spatial Signatures of Mouse Embryonic Endothelial Cells at Single-Cell Resolution.

BACKGROUND: Mature endothelial cells (ECs) are heterogeneous, with subtypes defined by tissue origin and position within the vascular bed (ie, artery, capillary, vein, and lymphatic). How this heterogeneity is established during the development of the vascular system, especially arteriovenous specification of ECs, remains incompletely characterized. METHODS: We used droplet-based single-cell RNA sequencing and multiplexed error-robust fluorescence in situ hybridization to define EC and EC progenitor subtypes from E9.5, E12.5, and E15.5 mouse embryos. We used trajectory inference to analyze the specification of arterial ECs (aECs) and venous ECs (vECs) from EC progenitors. Network analysis identified candidate transcriptional regulators of arteriovenous differentiation, which we tested by CRISPR (clustered regularly interspaced short palindromic repeats) loss of function in human-induced pluripotent stem cells undergoing directed differentiation to aECs or vECs (human-induced pluripotent stem cell-aECs or human-induced pluripotent stem cell-vECs). RESULTS: From the single-cell transcriptomes of 7682 E9.5 to E15.5 ECs, we identified 19 EC subtypes, including Etv2+Bnip3+ EC progenitors. Spatial transcriptomic analysis of 15 448 ECs provided orthogonal validation of these EC subtypes and established their spatial distribution. Most embryonic ECs were grouped by their vascular-bed types, while ECs from the brain, heart, liver, and lung were grouped by their tissue origins. Arterial (Eln, Dkk2, Vegfc, and Egfl8), venous (Fam174b and Clec14a), and capillary (Kcne3) marker genes were identified. Compared with aECs, embryonic vECs and capillary ECs shared fewer markers than their adult counterparts. Early capillary ECs with venous characteristics functioned as a branch point for differentiation of aEC and vEC lineages. CONCLUSIONS: Our results provide a spatiotemporal map of embryonic EC heterogeneity at single-cell resolution and demonstrate that the diversity of ECs in the embryo arises from both tissue origin and vascular-bed position. Developing aECs and vECs share common venous-featured capillary precursors and are regulated by distinct transcriptional regulatory networks.

A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network.

BACKGROUND: The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice. METHODS: FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription. RESULTS: FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency. CONCLUSIONS: Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

In vivo proximity proteomics uncovers palmdelphin (PALMD) as a Z-line-associated mitigator of isoproterenol-induced cardiac injury.

Z-lines are core ultrastructural organizers of cardiomyocytes that modulate many facets of cardiac pathogenesis. Yet a comprehensive proteomic atlas of Z-line-associated components remain incomplete. Here, we established an adeno-associated virus (AAV)-delivered, cardiomyocyte-specific, proximity-labeling approach to characterize the Z-line proteome in vivo. We found palmdelphin (PALMD) as a novel Z-line-associated protein in both adult murine cardiomyocytes and human pluripotent stem cell-derived cardiomyocytes. Germline and cardiomyocyte-specific palmd knockout mice were grossly normal at baseline but exhibited compromised cardiac hypertrophy and aggravated cardiac injury upon long-term isoproterenol treatment. By contrast, cardiomyocyte-specific PALMD overexpression was sufficient to mitigate isoproterenol-induced cardiac injury. PALMD ablation perturbed transverse tubules (T-tubules) and their association with sarcoplasmic reticulum, which formed the Z-line-associated junctional membrane complex (JMC) essential for calcium handling and cardiac function. These phenotypes were associated with disrupted localization of T-tubule markers caveolin-3 (CAV3) and junctophilin-2 (JPH2) and the reduction of nexilin (NEXN) protein, a crucial Z-line-associated protein that is essential for both Z-line and JMC structures and functions. PALMD was found to interact with NEXN and enhance its protein stability while the Nexn mRNA level was not affected. Together, this study discovered PALMD as a potential target for myocardial protection and highlighted in vivo proximity proteomics as a powerful approach to nominate novel players regulating cardiac pathogenesis. Highlights: In vivo proximity proteomics uncover novel Z-line components that are undetected in in vitro proximity proteomics in cardiomyocytes.PALMD is a novel Z-line-associated protein that is dispensable for baseline cardiomyocyte function in vivo.PALMD mitigates cardiac dysfunction and myocardial injury after repeated isoproterenol insults.PALMD stabilizes NEXN, an essential Z-line-associated regulator of the junctional membrane complex and cardiac systolic function.

Reduced Mitochondrial Protein Translation Promotes Cardiomyocyte Proliferation and Heart Regeneration.

BACKGROUND: The importance of mitochondria in normal heart function are well recognized and recent studies have implicated changes in mitochondrial metabolism with some forms of heart disease. Previous studies demonstrated that knockdown of the mitochondrial ribosomal protein S5 (MRPS5) by small interfering RNA (siRNA) inhibits mitochondrial translation and thereby causes a mitonuclear protein imbalance. Therefore, we decided to examine the effects of MRPS5 loss and the role of these processes on cardiomyocyte proliferation. METHODS: We deleted a single allele of MRPS5 in mice and used left anterior descending coronary artery ligation surgery to induce myocardial damage in these animals. We examined cardiomyocyte proliferation and cardiac regeneration both in vivo and in vitro. Doxycycline treatment was used to inhibit protein translation. Heart function in mice was assessed by echocardiography. Quantitative real-time polymerase chain reaction and RNA sequencing were used to assess changes in transcription and chromatin immunoprecipitation (ChIP) and BioChIP were used to assess chromatin effects. Protein levels were assessed by Western blotting and cell proliferation or death by histology and terminal deoxynucleotidyl transferase dUTP nick-end labeling (TUNEL) assays. Adeno-associated virus was used to overexpress genes. The luciferase reporter assay was used to assess promoter activity. Mitochondrial oxygen consumption rate, ATP levels, and reactive oxygen species were also analyzed. RESULTS: We determined that deletion of a single allele of MRPS5 in mice results in elevated cardiomyocyte proliferation and cardiac regeneration; this observation correlates with improved cardiac function after induction of myocardial infarction. We identified ATF4 (activating transcription factor 4) as a key regulator of the mitochondrial stress response in cardiomyocytes from Mrps5+/- mice; furthermore, ATF4 regulates Knl1 (kinetochore scaffold 1) leading to an increase in cytokinesis during cardiomyocyte proliferation. The increased cardiomyocyte proliferation observed in Mrps5+/- mice was attenuated when one allele of Atf4 was deleted genetically (Mrps5+/-/Atf4+/-), resulting in the loss in the capacity for cardiac regeneration. Either MRPS5 inhibition (or as we also demonstrate, doxycycline treatment) activate a conserved regulatory mechanism that increases the proliferation of human induced pluripotent stem cell-derived cardiomyocytes. CONCLUSIONS: These data highlight a critical role for MRPS5/ATF4 in cardiomyocytes and an exciting new avenue of study for therapies to treat myocardial injury.

Spatiotemporal cell junction assembly in human iPSC-CM models of arrhythmogenic cardiomyopathy.

Arrhythmogenic cardiomyopathy (ACM) is an inherited cardiac disorder that causes life-threatening arrhythmias and myocardial dysfunction. Pathogenic variants in Plakophilin-2 (PKP2), a desmosome component within specialized cardiac cell junctions, cause the majority of ACM cases. However, the molecular mechanisms by which PKP2 variants induce disease phenotypes remain unclear. Here we built bioengineered platforms using genetically modified human induced pluripotent stem cell-derived cardiomyocytes to model the early spatiotemporal process of cardiomyocyte junction assembly in vitro. Heterozygosity for truncating variant PKP2R413X reduced Wnt/ß-catenin signaling, impaired myofibrillogenesis, delayed mechanical coupling, and reduced calcium wave velocity in engineered tissues. These abnormalities were ameliorated by SB216763, which activated Wnt/ß-catenin signaling, improved cytoskeletal organization, restored cell junction integrity in cell pairs, and improved calcium wave velocity in engineered tissues. Together, these findings highlight the therapeutic potential of modulating Wnt/ß-catenin signaling in a human model of ACM.

Fibre-infused gel scaffolds guide cardiomyocyte alignment in 3D-printed ventricles.

Hydrogels are attractive materials for tissue engineering, but efforts to date have shown limited ability to produce the microstructural features necessary to promote cellular self-organization into hierarchical three-dimensional (3D) organ models. Here we develop a hydrogel ink containing prefabricated gelatin fibres to print 3D organ-level scaffolds that recapitulate the intra- and intercellular organization of the heart. The addition of prefabricated gelatin fibres to hydrogels enables the tailoring of the ink rheology, allowing for a controlled sol-gel transition to achieve precise printing of free-standing 3D structures without additional supporting materials. Shear-induced alignment of fibres during ink extrusion provides microscale geometric cues that promote the self-organization of cultured human cardiomyocytes into anisotropic muscular tissues in vitro. The resulting 3D-printed ventricle in vitro model exhibited biomimetic anisotropic electrophysiological and contractile properties.

Tbx5 maintains atrial identity by regulating an atrial enhancer network.

Understanding how the atrial and ventricular chambers of the heart maintain their distinct identity is a prerequisite for treating chamber-specific diseases. Here, we selectively inactivated the transcription factor Tbx5 in the atrial working myocardium of the neonatal mouse heart to show that it is required to maintain atrial identity. Atrial Tbx5 inactivation downregulated highly chamber specific genes such as Myl7 and Nppa , and conversely, increased the expression of ventricular identity genes including Myl2 . Using combined single nucleus transcriptome and open chromatin profiling, we assessed genomic accessibility changes underlying the altered atrial identity expression program, identifying 1846 genomic loci with greater accessibility in control atrial cardiomyocytes compared to KO aCMs. 69% of the control-enriched ATAC regions were bound by TBX5, demonstrating a role for TBX5 in maintaining atrial genomic accessibility. These regions were associated with genes that had higher expression in control aCMs compared to KO aCMs, suggesting they act as TBX5-dependent enhancers. We tested this hypothesis by analyzing enhancer chromatin looping using HiChIP and found 510 chromatin loops that were sensitive to TBX5 dosage. Of the loops enriched in control aCMs, 73.7% contained anchors in control-enriched ATAC regions. Together, these data demonstrate a genomic role for TBX5 in maintaining the atrial gene expression program by binding to atrial enhancers and preserving tissue-specific chromatin architecture of atrial enhancers.

Dynamic changes in P300 enhancers and enhancer-promoter contacts control mouse cardiomyocyte maturation.

Cardiomyocyte differentiation continues throughout murine gestation and into the postnatal period, driven by temporally regulated expression changes in the transcriptome. The mechanisms that regulate these developmental changes remain incompletely defined. Here, we used cardiomyocyte-specific ChIP-seq of the activate enhancer marker P300 to identify 54,920 cardiomyocyte enhancers at seven stages of murine heart development. These data were matched to cardiomyocyte gene expression profiles at the same stages and to Hi-C and H3K27ac HiChIP chromatin conformation data at fetal, neonatal, and adult stages. Regions with dynamic P300 occupancy exhibited developmentally regulated enhancer activity, as measured by massively parallel reporter assays in cardiomyocytes in vivo, and identified key transcription factor-binding motifs. These dynamic enhancers interacted with temporal changes of the 3D genome architecture to specify developmentally regulated cardiomyocyte gene expressions. Our work provides a 3D genome-mediated enhancer activity landscape of murine cardiomyocyte development.

Genetic modifiers modulate phenotypic expression of tafazzin deficiency in a mouse model of Barth syndrome.

Barth syndrome is an X-linked disorder caused by loss-of-function mutations in Tafazzin (TAZ), an acyltransferase that catalyzes remodeling of cardiolipin, a signature phospholipid of the inner mitochondrial membrane. Patients develop cardiac and skeletal muscle weakness, growth delay and neutropenia, although phenotypic expression varies considerably between patients. Taz knockout mice recapitulate many of the hallmark features of the disease. We used mouse genetics to test the hypothesis that genetic modifiers alter the phenotypic manifestations of Taz inactivation. We crossed TazKO/X females in the C57BL6/J inbred strain to males from eight inbred strains and evaluated the phenotypes of first-generation (F1) TazKO/Y progeny, compared to TazWT/Y littermates. We observed that genetic background strongly impacted phenotypic expression. C57BL6/J and CAST/EiJ[F1] TazKO/Y mice developed severe cardiomyopathy, whereas A/J[F1] TazKO/Y mice had normal heart function. C57BL6/J and WSB/EiJ[F1] TazKO/Y mice had severely reduced treadmill endurance, whereas endurance was normal in A/J[F1] and CAST/EiJ[F1] TazKO/Y mice. In all genetic backgrounds, cardiolipin showed similar abnormalities in knockout mice, and transcriptomic and metabolomic investigations identified signatures of mitochondrial uncoupling and activation of the integrated stress response. TazKO/Y cardiac mitochondria were small, clustered and had reduced cristae density in knockouts in severely affected genetic backgrounds but were relatively preserved in the permissive A/J[F1] strain. Gene expression and mitophagy measurements were consistent with reduced mitophagy in knockout mice in genetic backgrounds intolerant of Taz mutation. Our data demonstrate that genetic modifiers powerfully modulate phenotypic expression of Taz loss-of-function and act downstream of cardiolipin, possibly by altering mitochondrial quality control.

In Vivo Dissection of Chamber-Selective Enhancers Reveals Estrogen-Related Receptor as a Regulator of Ventricular Cardiomyocyte Identity.

BACKGROUND: Cardiac chamber-selective transcriptional programs underpin the structural and functional differences between atrial and ventricular cardiomyocytes (aCMs and vCMs). The mechanisms responsible for these chamber-selective transcriptional programs remain largely undefined. METHODS: We nominated candidate chamber-selective enhancers (CSEs) by determining the genome-wide occupancy of 7 key cardiac transcription factors (GATA4, MEF2A, MEF2C, NKX2-5, SRF, TBX5, TEAD1) and transcriptional coactivator P300 in atria and ventricles. Candidate enhancers were tested using an adeno-associated virus-mediated massively parallel reporter assay. Chromatin features of CSEs were evaluated by performing assay of transposase accessible chromatin sequencing and acetylation of histone H3 at lysine 27-HiChIP on aCMs and vCMs. CSE sequence requirements were determined by systematic tiling mutagenesis of 29 CSEs at 5 bp resolution. Estrogen-related receptor (ERR) function in cardiomyocytes was evaluated by Cre-loxP-mediated inactivation of ERRα and ERRγ in cardiomyocytes. RESULTS: We identified 134 066 and 97 506 regions reproducibly occupied by at least 1 transcription factor or P300, in atria or ventricles, respectively. Enhancer activities of 2639 regions bound by transcription factors or P300 were tested in aCMs and vCMs by adeno-associated virus-mediated massively parallel reporter assay. This identified 1092 active enhancers in aCMs or vCMs. Several overlapped loci associated with cardiovascular disease through genome-wide association studies, and 229 exhibited chamber-selective activity in aCMs or vCMs. Many CSEs exhibited differential chromatin accessibility between aCMs and vCMs, and CSEs were enriched for aCM- or vCM-selective acetylation of histone H3 at lysine 27-anchored loops. Tiling mutagenesis of 29 CSEs identified the binding motif of ERRα/γ as important for ventricular enhancer activity. The requirement of ERRα/γ to activate ventricular CSEs and promote vCM identity was confirmed by loss of the vCM gene profile in ERRα/γ knockout vCMs. CONCLUSIONS: We identified 229 CSEs that could be useful research tools or direct therapeutic gene expression. We showed that chamber-selective multi-transcription factor, P300 occupancy, open chromatin, and chromatin looping are predictive features of CSEs. We found that ERRα/γ are essential for maintenance of ventricular identity. Finally, our gene expression, epigenetic, 3-dimensional genome, and enhancer activity atlas provide key resources for future studies of chamber-selective gene regulation.

Ryanodine receptor 2 (RYR2) dysfunction activates the unfolded protein response and perturbs cardiomyocyte maturation.

AIMS: Calcium-handling capacity is a major gauge of cardiomyocyte maturity. Ryanodine receptor 2 (RYR2) is the pre-dominant calcium channel that releases calcium from the sarcoplasmic reticulum/endoplasmic reticulum (SR/ER) to activate cardiomyocyte contraction. Although RYR2 was previously implied as a key regulator of cardiomyocyte maturation, the mechanisms remain unclear. The aim of this study is to solve this problem. METHODS AND RESULTS: We performed Cas9/AAV9-mediated somatic mutagenesis to knockout RYR2 specifically in cardiomyocytes in mice. We conducted a genetic mosaic analysis to dissect the cell-autonomous function of RYR2 during cardiomyocyte maturation. We found that RYR2 depletion triggered ultrastructural and transcriptomic defects relevant to cardiomyocyte maturation. These phenotypes were associated with the drastic activation of ER stress pathways. The ER stress alleviator tauroursodeoxycholic acid partially rescued the defects in RYR2-depleted cardiomyocytes. Overexpression of ATF4, a key ER stress transcription factor, recapitulated defects in RYR2-depleted cells. Integrative analysis of RNA-Seq and bioChIP-Seq data revealed that protein biosynthesis-related genes are the major direct downstream targets of ATF4. CONCLUSION: RYR2-regulated ER homeostasis is essential for cardiomyocyte maturation. Severe ER stress perturbs cardiomyocyte maturation primarily through ATF4 activation. The major downstream effector genes of ATF4 are related to protein biosynthesis.

Tbx5 maintains atrial identity in post-natal cardiomyocytes by regulating an atrial-specific enhancer network.

Understanding how the atrial and ventricular heart chambers maintain distinct identities is a prerequisite for treating chamber-specific diseases. Here, we selectively knocked out (KO) the transcription factor Tbx5 in the atrial working myocardium to evaluate its requirement for atrial identity. Atrial Tbx5 inactivation downregulated atrial cardiomyocyte (aCM) selective gene expression. Using concurrent single nucleus transcriptome and open chromatin profiling, genomic accessibility differences were identified between control and Tbx5 KO aCMs, revealing that 69% of the control-enriched ATAC regions were bound by TBX5. Genes associated with these regions were downregulated in KO aCMs, suggesting they function as TBX5-dependent enhancers. Comparing enhancer chromatin looping using H3K27ac HiChIP identified 510 chromatin loops sensitive to TBX5 dosage, and 74.8% of control-enriched loops contained anchors in control-enriched ATAC regions. Together, these data demonstrate TBX5 maintains the atrial gene expression program by binding to and preserving the tissue-specific chromatin architecture of atrial enhancers.

RBPMS2 Is a Myocardial-Enriched Splicing Regulator Required for Cardiac Function.

BACKGROUND: RBPs (RNA-binding proteins) perform indispensable functions in the post-transcriptional regulation of gene expression. Numerous RBPs have been implicated in cardiac development or physiology based on gene knockout studies and the identification of pathogenic RBP gene mutations in monogenic heart disorders. The discovery and characterization of additional RBPs performing indispensable functions in the heart will advance basic and translational cardiovascular research. METHODS: We performed a differential expression screen in zebrafish embryos to identify genes enriched in nkx2.5-positive cardiomyocytes or cardiopharyngeal progenitors compared to nkx2.5-negative cells from the same embryos. We investigated the myocardial-enriched gene RNA-binding protein with multiple splicing (variants) 2 [RBPMS2)] by generating and characterizing rbpms2 knockout zebrafish and human cardiomyocytes derived from RBPMS2-deficient induced pluripotent stem cells. RESULTS: We identified 1848 genes enriched in the nkx2.5-positive population. Among the most highly enriched genes, most with well-established functions in the heart, we discovered the ohnologs rbpms2a and rbpms2b, which encode an evolutionarily conserved RBP. Rbpms2 localizes selectively to cardiomyocytes during zebrafish heart development and strong cardiomyocyte expression persists into adulthood. Rbpms2-deficient embryos suffer from early cardiac dysfunction characterized by reduced ejection fraction. The functional deficit is accompanied by myofibril disarray, altered calcium handling, and differential alternative splicing events in mutant cardiomyocytes. These phenotypes are also observed in RBPMS2-deficient human cardiomyocytes, indicative of conserved molecular and cellular function. RNA-sequencing and comparative analysis of genes mis-spliced in RBPMS2-deficient zebrafish and human cardiomyocytes uncovered a conserved network of 29 ortholog pairs that require RBPMS2 for alternative splicing regulation, including RBFOX2, SLC8A1, and MYBPC3. CONCLUSIONS: Our study identifies RBPMS2 as a conserved regulator of alternative splicing, myofibrillar organization, and calcium handling in zebrafish and human cardiomyocytes.